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Isolated Millimeter Flares of Cyg X-3 
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Abstract 

Cygnus X-3 (Cyg X-3) is a well-known microquasar with rclativistic jets. Cyg X-3 is especially famous 
for its giant radio outbursts, which have been observed once every few years since their first discovery. Each 
giant outburst presumably consists of a series of short-duration flares. The physical parameters of the flares 
in the giant outbursts are difficult to derive because the successive flares overlap. Here, we report isolated 
flares in the quiescent phase of Cyg X-3, as observed at 23, 43, and 86 GHz with the 45-m radio telescope 
at Nobeyama Radio Observatory. The observed flares have small amplitude (0.5-2 Jy) and short duration 
(1-2 h). The millimeter fluxes rapidly increase and then exponentially decay. The lifetime of the decay 
is shorter at higher frequency. The radio spectrum of Cyg X-3 during the flares is flat or inverted around 
the peak flux density. After that, the spectrum gradually becomes steeper. The observed characteristics 
are consistent with those of adiabatic expanding plasma. The brightness temperature of the plasma at the 
peak is estimated to be Tb ^ 1 x 10^^ K. The magnetic field in the plasma is calculated to be 0.2 < ^ 30 
G. 
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1. Introduction 

Cygnus X-3 (Cyg X-3) is a well-known microquasar that 
contains a Wolf-Rayet star and a compact object, the na- 
ture and mass of which are still not known (Lommen et al. 
2005;Shrader et al 2010). Giant radio outbursts from Cyg 
X-3 have been observed once every few years since their 
first discovery. Very long baseline interferometry observa- 
tions of these outbursts revealed that Cyg X-3 sometimes 
has a bipolar relativistic jet (Marti et al. 2001, Miller- 
Jones et al. 2004, Tudose et al. 2007) and sometimes 
has a one-sided jet (Mioduszewski et al. 2001). These 
relativistic jets are presumably the remnants of ejection 
from the flares. In addition, Cyg X-3 was first detected 
in the GeV 7-ray regime in 2009 by the Fermi (Abdo 
et al. 2009) and AGILE (Tavani et al. 2009) satellites. 
Williams et al. reported that a radio flare was followed by 
a 7-ray flare (Williams et al. 2011). The close relation- 
ship between radio and 7-ray flares is also supported by 
new observations (Kotani et al. 2011, Bulgarelli et al. 
2011). The centimeter wave emission from Cyg X-3 has 
been monitored by several groups, even during the quies- 
cent phase (e.g. Trushkin, Nizhelskij, & Bursov 2008). 
The Green Bank interferometer at 5 GHz and the Ryle 
telescope at 15 GHz have revealed the detailed behavior 



of radio emission from Cyg X-3 with a high time reso- 
lution (e.g. Ogley et al. 2001). The centimeter radio 
emission changes gradually over a time scale of 1 h or 
less. In contrast, the time scale of the change during out- 
bursts is shorter at higher frequency. In particular, the 
time scale for millimeter waves is much shorter than that 
for centimeter waves (e.g. Tsuboi et al. 2008). 

The giant outbursts of Cyg X-3 are presumed to consist 
of a series of short-duration flares. In previous research, a 
millimeter wave monitor with sufficient time coverage and 
high time resolution found a new transient phenomenon of 
Cyg X-3: a flux density change with an e- folding rise time 
of 3.6 min or less in the millimeter light curve (Tsuboi 
et al. 2010). However, such flares in the giant outburst 
are too crowded together to discriminate the physical pa- 
rameters of a single flare, even in the millimeter wave re- 
gion. Ascertaining these parameters is vital for exploring 
the mechanism of the flares. 

Accordingly, to obtain light curves that do not suf- 
fer from the interference of successive flares, observa- 
tions were performed in the millimeter wave region at a 
high time resolution during quiescent phase of Cyg X-3 
by using the 45-m radio telescope at Nobeyama Radio 
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Observatory^ (NR045). This telescope has wide fre- 
quency coverage from 20 to 100 GHz. Studying the be- 
havior of flares at different frequencies also provided im- 
portant information for exploring the mechanism of the 
flares. 

2. Observations 

We monitored the flux density of Cyg X-3 in the quies- 
cent phase at a high time resolution (4-10 min) by using 
NR045 at 23, 43, and 86 GHz for 20 days during the 
period from 2006 to 2010. 

To conduct the observations, a cooled high electron mo- 
bility transistor receiver with a dual circular polarization 
feed was used at 23.0 GHz. However, observations at 
23 GHz were performed only in 2006 because changing 
the frequency between 23 GHz and 43 and 86 GHz re- 
quired a dead time of 5 min, which we judged to be a 
long time when compared with the time scale of Cyg X- 
3 phenomena. Superconductor-insulator-superconductor 
receivers with orthogonal linear polarization feeds were 
used to make observations at 43.0 and 86.0 GHz. The 
frequency bandwidth of these receivers is 500 MHz. The 
system noise temperatures during the observations, in- 
cluding atmospheric effects and antenna ohmic loss, were 
80-120 K at 23 GHz, 120-200 K at 43 GHz, and 250-350 
K at 86 GHz. The full width at half- maximum values of 
the telescope beams were 77" at 23 GHz, 39" at 43 GHz, 
and 19" at 86 GHz. To remove any atmospheric effects, 
the telescope beam was alternated at 15 Hz between the 
position of the source and a region of blank sky 6' from 
the source in the azimuth direction by a beam switch. The 
telescope was position-switched between an on-source po- 
sition at A / cos El = 0' and a region of blank sky at 
AAz/ cos El ~ 3', where the duration at each position was 
8 s/cycle. The antenna temperatures were calibrated by 
a chopper wheel method. The primary flux calibrator for 
the conversion from antenna temperature to flux density 
was a proto-planetary nebula, NGC 7027, whose flux den- 
sity values are given as 5.5 Jy at 22 GHz, 5.0 Jy at 43 
GHz, and 4.6 Jy at 86 GHz (Ott et al. 1994). Telescope 
pointing was checked and corrected in every observation 
procedure by viewing NGC 7027 in cross-scan mode. The 
pointing accuracy was better than 3" root mean square 
during the observations. 

The uncertainty in the flux density of Cyg X-3 depends 
on the weather conditions. However, the sensitivity of 
the telescopes is not a principal factor in the uncertainty. 
Because the primary flux-scale calibrator, NGC 7027, is 
close to Cyg X-3 in the celestial sphere and the differ- 
ence in atmospheric attenuation between these sources 
is calibrated by the chopper wheel method, there is no 
significant effect on the flux density data. The typical 
systematic uncertainty for NR045 is 10%; however, the 
relative uncertainties in the flux density on any particular 
day should be much less than this value because the ob- 



servations were performed under fine weather conditions. 
3. Results 

From the 20 days of monitoring during the quiescent 
phase of Cyg X-3, we found two isolated flares, at MJD 
= 53809 (15 March 2006 UT) and at MJD = 54972 (May 
20 2009 UT). We observed the whole episode of these flares 
without contamination from successive flares. There was 
no isolated flare with statistical significance on other days. 
The complete monitoring data will be published in a fu- 
ture paper. 

3.1. Flare at MJD 53809 

Figure la shows the light curves of the flare at 
MJD53809 at 23, 43, and 86 GHz. The MJD =53809 flare 
was found in the post-bursting phase of the giant outburst 
in spring 2006. The flux densities rapidly increased to a 
peak around MJD = 53809.78. Unfortunately, because 
calibration was performed near the peak, there is a break 
in the data. The maximum flux densities were at least 
= 1.17 ±0.01 Jy at 43 GHz and 5"^ = 1.71 ±0.05 Jy 
at 86 GHz, where the error includes only statistical error. 
The time scales of the e- folding rise, te, were below tp ^ 7 
min at 43 GHz and ^6 min at 86 GHz. In an earlier 
paper (Tsuboi et al. 2010), we reported an abrupt large- 
amplitude change in flux density with an e-folding rise 
time of 4 min or less in light curves at 100 GHz (Tsuboi 
et al. 2010). The time scales shown in Figure la are simi- 
lar to the previous value. In contrast, the flux density ap- 
pears to change more gradually at 23 GHz than at 43 GHz 
and 86 GHz, suggesting that the time scale becomes longer 
at lower frequency. However, the time lag between 43 and 
86 GHz is not clear. The absence of time lag is similar to 
that seen for light curves at 100 GHz (Tsuboi et al. 2010). 
After the intensity peak, the flux densities then exponen- 
tially decayed. The lifetimes (5,^ oc cxp { — {t ~ t^) / t)) of 
the fluxes were r = 0.167to °2i days at 23 GHz (dotted 
line in Figure la), r = O.OdStoMil days at 43 GHz (bro- 
ken line in Figure la), and t = 0.033to:oo5 days at 86 GHz 
(solid line in Figure la). Here, the initial time of the flare, 
^0, is assumed to be the same as the time at the peak flux 
density, to = 53809.78 in MJD. The hfetime of the flux 
was much shorter at higher frequency. In addition, there 
might be a plateau around MJD = 53809.83-53809.86 in 
the flux density at 86 GHz. A plateau was not clearly seen 
in the flux densities at 23 and 43 GHz. 

Figure 2a shows the radio spectral index of Cyg X-3 at 
MJD = 53809. The values are derived from the inclination 
between the flux densities at 43 and 86 GHz. The radio 
spectrum of Cyg X-3 is inverted, a ~ —0.5(5*1/ oc i^~"), 
around the peak flux density. After the peak flux density, 
the spectral index rapidly becomes steeper; however, dur- 
ing the exponential decay of the flux density, the spectral 
index is around a ^ 0.5. 
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3.2. Flare at MJD = 54972 

Figure lb shows the hght curves of the flare at MJD 
= 54972 at 43 and 86 GHz. MJD = 54972 corresponded 
to the quiescent phase of Cyg X-3 in spring 2009. This 
flare has been briefly reported (Tsuboi et al. 2009). The 
flux densities rapidly increased to a peak around MJD = 
54972.88. The maximum flux densities are = 0.43±0.02 
Jy at 43 GHz and = 0.39 ± 0.02 Jy at 86 GHz. The 
time scales of the e- folding rise are ig ^ 11 min at 43 GHz 
and ie ^ 7 min at 86 GHz, slightly longer than those in 
Figure la. The flux densities exponentially decayed after 
the intensity peaked. The lifetimes of the fluxes, assuming 
to = 54972.88, are t = 0.077l°:oo3 days at 43 GHz (broken 
line in Figure lb), and t = O.OSS+^^J^ days at 86 GHz 
(solid line in Figure lb). These lifetimes are consistent 
with those in Figure la. There might also be a plateau 
in the flux density at 86 GHz around MJD = 54972.91- 
54972.95. In contrast to the first flare, a similar feature is 
also observed at 43 GHz. 

Figure 2b shows the radio spectral index of Cyg X-3 
at MJD = 54972. The radio spectrum of Cyg X-3 is flat 
around the peak flux density, but rapidly become steeper 
to a ~ 0.5 during the exponential decay of the flux density. 
This is consistent with the spectral index in Figure 2a. 
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Fig. 1. a Light curves of Cyg X-3 at MJD = 53809 ob- 
served with the NR045 telescope. The frequency bands are 
23 (cross), 43 (open circle), and 86 (filled circle) GHz. The 
maximum flux densities occur around MJD = 53809.78. They 
are = l.f7±0.0f Jy at 43 GHz and S„ = 1.71 ±0.05 Jy 
at 86 GHz. After the intensity peak, the flux densities expo- 
nentially decayed. The lifetimes (5v oc cxp (— (i — to)/T)) are 
r = 0.167^0 021 '^^y^ 23 GHz (dotted line), t = 0.098 
days at 43 GHz (broken line), and r = 0.033lo "os 
86 GHz (solid line). Here, the initial time of the flare, to, 
is assumed to be the time at the peak, to = 53809.78. b 
Light curves of Cyg X-3 at MJD = 54972 observed with the 
NR045 telescope. The frequency bands are 43 (open circle) 
and 86 (filled circle) GHz. The maximum flux densities occur 
around MJD = 54972.88. They are 5^ = 0.43 ± 0.02 Jy at 
43 GHz and = 0.39 ± 0.02 Jy at 86 GHz. The lifetimes. 
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Fig. 2. a Spectral index between 43 and 86 GHz (filled cir- 
cle) of Cyg X-3 at MJD = 53809 observed with the NR045 
telescope. The flux density at 86 GHz (open circle) is also 
shown for comparison. The radio spectrum of Cyg X-3 is 
inverted, a ~ —0.5(5,^ oc around the peak flux den- 

sity. The spectral index was around a ~ 0.5 during the ex- 
ponential decay of the flux density, b Spectral index be- 
tween 43 and 86 GHz (filled circle) of Cyg X-3 at MJD 
= 54972 observed with the NR045 telescope. The flux den- 
sity at 86 GHz (open circle) is also shown for compari- 
son. The radio spectrum of Cyg X-3 has a flat spectrum 
around the peak flux density. The spectral index was at 
around a 0.5 during the exponential decay of the flux. 
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4. Discussion 

The time lag between 43 and 86 GHz was not observed 
in the rising phase of the flares; the two flares had flat 
or inverted spectra around the maximum intensity. The 
flux densities of both flares exponentially decayed over 
a similar lifetime and with a similar spectral index after 
the peak. These observed characteristics are consistent 
with those of adiabatic expanding plasma; that is, an op- 
tically thick plasma blob is ejected. This blob rapidly 
expands and its flux density increases rapidly, indepen- 
dent of frequency. The optical thickness decreases with 
the expansion and the plasma blob becomes optically thin 
at a higher frequency. The magnetic field producing the 
synchrotron emission decays. The emission at higher fre- 
quency decays more quickly. 

The radius of the emitting region at the maximum in- 
tensity, r, can be estimated from r < cie/2, and the upper 
limit of the source size is r < 0.4 AU. The observation fre- 
quency is 1/ = 86 GHz and the flux density is 5',^ = 1.7 Jy 
for the flare at MJD — 53809. The brightness temperature 
of the emitting region, Tb, is given by 



(1) 



where d is the distance of Cyg X-3 and fc^ is Boltzmann 
conatatnt (Ogley et al. 2001). This brightness tempera- 
ture is larger than Tb > 1.5 x 10^^ K, assuming that Cyg 
X-3 is located at a distance of d = 9 kpc. 

The observed spectrum of the flares around the peak is 
flat or inverted between i/ = 43 and 86 GHz. If the emission 
is completely optically thick in the frequency range, the 
spectral index should be a ^ —2.5 for relativistic electrons 
with a power energy spectrum, N{E) oc E~'^. The ob- 
served spectral indexes indicate that emission in this fre- 
quency range is close to the turnover frequency, at which 
the optical thickness becomes r ~ 1 and is given by 



"^14 w) — t7:jt^ 



(2) 



(seeEq. 6.38 in Pacholczyk 1970). 

We assume that the electron energy index is 7 = 2, which 
is derived from the spectral index of a ^ 0.5 during the ex- 
ponential decay. Hence, with 7 = 2 the turnover frequency 
becomes 

I'^lGHz] -42 X iJj_[G]i (3) 

Here, we use the coefficient Ci4(2) — 3.42 x lO'^". On the 
other hand, the turnover frequency in the flares is esti- 
mated to be lyg ^ 80 GHz. owing to the observed flares 
showing a flat or inverted spectrum between ;/ = 43 and 
86 GHz. The magnetic field strength at the peak is calcu- 
lated to be ~ 30 G and the cooling time of the electrons, 
tc, is estimated to be 



(4) 



where C12 = 1.6 x 10^ (sec Eq. 7.7 in Pacholczyk 1970). 
If the magnetic field is H± < 30 G, the cooling time is 



tc ^ 14 h, which is consistent with no change occurring in 
the spectral index during the observed exponential decay. 
The exponential decay is assumed to be caused by the 
decrease in the strength of the magnetic field. 

Assuming equipartition between the magnetic field and 
the electrons, the minimum energy of the source to ac- 
count for the observed synchrotron emission is approxi- 
mated by 



En 



3 X lO-'^ 1 
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9/7 



GHz 



VmJy, 



2.1 X 10-*" erg 



where Cp/eo is the ratio between the energy in protons 
and that in electrons; this ratio is assumed to be zero 
hereafter (Fender 2003). The minimum power given to 
the electrons is found by dividing the minimum energy by 
the time scale of the e- folding rise to- 

P.„i„- ^=5.1x10=^7 erg s-i. (7) 

The power is comparable with those estimated from flux 
densities at low frequencies, whose rise time scales are 
longer. 

From the assumption of equipartition, the magnetic 
field strength, H^q, is 
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Note that the equipartition magnetic field strength, H^q ^ 
10 G, is of the same order as the strength of ~ 30 G, 
obtained from the turnover frequency. 

The ratio of the inverse-Compton scattering loss, 

, and the synchrotron radiation loss, i for 

the relativistic electrons is given by 



R- 



'^rad 
^mag 



Si 



(9) 



where Urad and Umag are the radiation energy density and 
the magnetic energy density, respectively (see Eq. 5.60 in 
Pacholczyk 1970, Tsuboi et al. 2010). If the inverse- 
Compton scattering loss dominates the synchrotron radi- 
ation loss, the lifetime of the relativistic electron would be 
so short that the emission would be invisible. Therefore, 
we set the ratio i? < 1 by assuming that the spectral in- 
dex is flat until 1/ — 100 GHz. In this case, the magnetic 
field is greater than H > 0.2 G and thus, the magnetic 
field strength of Cyg X-3 is estimated to be in the range 
0.2<i? <30 G. 

5. Conclusions 

We performed millimeter wave monitor observations of 
Cyg X-3 in quiescent phases between giant outbursts at 
23, 43, and 86 GHz by using the NR045 radio telescope. 
We found two isolated flares with amplitude 0.5-2 Jy, for 
which the light curves did not suffer from successive flares. 
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The e-folding rise times were as short as 6-7 min at 86 
GHz. The spectra were fiat or slightly inverted around 
the maximum intensity. The flux densities of both flares 
exponentially decayed over a similar lifetime, r = 0.03 days 
at 86 GHz, and with a similar spectral index, a ^ 0.5. The 
observed characteristics are consistent with those of adia- 
batic expanding plasma. The observed e-folding rise times 
showed that the upper limit radius of the ejecting plasma 
was smaller than r < 0.4 AU. The brightness temperature 
of the ejecting plasma was estimated to be Tb ^ 1 x 10^^ 
K. The magnetic field in the plasma was derived to be 
0.2<i? <30 G. 

The authors thank the members of the 45-m radio tele- 
scope group at Nobeyama Radio Observatory for sup- 
port in the observations. This research was supported by 
KAKENHI 19047002 (TK) and by KAKENHI 23340075 
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